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Abstract 

The effect of magnetic field on the nonlinear growth rate of Rayleigh - Taylor instability 
induced two fluid interfacial structures has been investigated. The magnetic field is assumed to 
be parallel to the plane of the two fluid interface and acts in a direction perpendicular to the wave 
vector. If magnetic field is restricted only to either side of the interface the growth rate may be 
depressed (may almost disappear) or be enhanced depending on whether the magnetic pressure 
on the interface opposes the instability driving pressure difference g{ph — Pi)y or acts in the 
same direction. If magnetic field is present on both sides of the two fiuid interface, stabilization 
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may also take place in the sense that the surface of separation undulates periodically when the 
force due to magnetic pressure on two sides are such as to act in opposite direction. This result 
differs from the classical linear theory result which predicts that the magnetic field parallel to 
the surface has no influence on the growth rate when the wave vector is perpendicular to its 
direction. 
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I. INTRODUCTION 

Temporal development of nonlinear structures at the two fluid interface consequent to Rayleigh - 
Taylor (RT) or Richtmyer - Meshkov (RM) instability is of much current interest both from theoretical 
and experimental point of view. The structure is called a bubble if the lighter fluid pushes across the 
unperturbed interface into the heavier fluid and a spike if the opposite takes place. The importance 
of such instabilities arises in connection with a wide range of problems ranging from astrophysical 
phenomena such as Supernova remnant to Inertial Confinement Fusion (ICF). A core collapse Super 
Nova (SN) is driven by an externally powerful shock, and strong shocks are the breeding ground of 
hydrodynamic instabilities like RT and RM instabilities. During the shock transit phase, the RM 
instability is activated at each discontinuity in the density profile of the star at the 0-He and He-H 
interface. After shock transit, hydrodynamic mixing continues due to RT instability, as the denser 
layers are decelerated by lower density outer layer. 

In an ICF situation, ablation front of an imploding capsule is subject to the RT instability because 
dense core is compressed and accelerated by low density ablating plasma. RT instability enhances 
the perturbation initiated by laser induced target non uniformity and consequently the performance 
of ICF implosion may be seriously affected. The dynamics of the instability of the interface of two 
constant density non-conducting fluids and the associated nonlinear structure has been studied by 
several authors using an expression near the tip of the bubble or the spike up to second order 

in the transverse coordinate following Layzer's approach The fluids may also be ionized as in the 
astrophysical situation or may get ionized through laser irradiation in laboratory condition. Magnetic 
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field generated by ponderomotive force can exist l^"!^ in sucli conducting (ionized) fluids and have 
important influence on the growth or suppression of the instabilities. 

When the explosion of a la type supernova (SNIa) starts in a white dwarf as a laminar deflagration 
at the center of the star, RT instability begins to act [IHl.III] The burning velocity at these regimes 
can be described by fractal model of combustion. In white dwarf, magnetic field with strength 
upto 10® ~ 10^ G exist at the surface and the field near the center may be ~ 10 times greater. 
Rayleigh - Taylor instability arising during type la supernova explosion is associated with strong 
magnetic field. Since the magnetic field is dipolar type the fluid propagates parallel to the field lines 
(i.e., approximately along the direction of gravity) near the magnetic pole while the field lines are 
transverse to the direction of gravity at the magnetic equatorial region. Thus magnetic field effect 
on RT instability may have important roles to play whether the field lines are normal or parallel to 
the two fluid interface (i.e., along or perpendicular to the direction of gravity). 

The effect of magnetic field on Rayleigh - Taylor instability has been studied in detail previously 
by Chandrasekhar When the magnetic field is normal to the surface of separation of the two 
fluids, the RT instability is almost unaffected by the magnetic field when the wave number 'k' of 
the perturbation is small; but contrary to the purely hydrodynamic case the growth rate does not 
increase indefinitely with 'k' but tends to a saturation value as — )■ oo. Magnetic field parallel to 

(TO' 

the direction of impulsively generated acceleration is also shown to induce RM instability. This 
however happens for sufficiently intense magnetic field and also tends asymptotically to a saturation 
value. 

In case the magnetic field is parallel to the surface of separation it is found that according to linear 
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theory there exists no effect of the magnetic field on the instabihty if the latter is perpendicular 
to the wave vector k. Non vanishing effect of transverse magnetic field with k perpendicular to the 
zeroeth order magnetic field is however found to exist in linear theory when compressibility effect is 



taken into account ESI. The growth rate is found to be lowered both for continuously accelerated 



(RTI) and impulsively accelerated (RMI) two fluid interface LL^-L^ when k has component parallel 
to the magnetic field. The nature of the depression has close resemblance to that due to surface 
tension ^2] ^^d also has useful application in astrophysical context M'lSS] 

The present paper is addressed to the problem of the time development of the nonlinear interfacial 
structure caused by Rayleigh Taylor instability in presence of a magnetic field parallel to the surface 
of separation of the two fluids . The wave vector is assumed to lie in the same plane and perpendicular 
to the magnetic field. With such a geometry there is no effect of the magnetic field in the classical 
12] linear approximation. However, it is no longer the case when linearization restriction is lifted. 
This may be understood from the following consideration. 

In presence of magnetic field, there exists the magnetic pressure in addition to the usual hydro- 
dynamic pressure. As a result the RT instability driving pressure difference g{ph — Pi)j is changed by 
the inclusion of the magnetic pressure difference {l/2p){Bf^ — Bf) [the suffix h(l) correspond to the 
dynamical variable associated with the heavier (lighter) fluid]. This has the consequence that the 
growth rate may be enhanced or depressed according as the extra contribution is either positive or 
negative . Moreover, as we shall see there may also occur stabilization in the sense that the surface 
of separation executes periodic undulation resulting from time lag in the temporal variation of Bh 
and Bi. It is interesting to note that these are entirely nonlinear effects and disappear in the linear 



approximation. 

Section II deals with the basic MHD equation together with the geometry involved. The fluid is 
assumed inviscid and perfectly conducting and the fluid motion to be one of potential type motion. 
The investigation of the nonlinear aspect of the mushroom structure of the two fluid interface is 
facilitated by Bernoulli's equation - the first integral of the equation of motion of the magnetofluid 
obtained with the help of the magnetic field geometry. The kinematical and dynamical boundary 
conditions holding at the two fluid interface are set forth in section III. The set of equations describing 
the temporal development of the RT instability induced nonlinear structures at the interface are 
derived in section IV. As these equation are not amenable to solution in closed analytic form, the 
results are obtained by numerical methods followed by graphical results and are presented in section 
V. A summary of the results is given in section VI. 



Assume that the undisturbed surface is y = 0, the transverse coordinates being represented by x, z. 
The heavier fluid (density ph =constant) occupies the region |/ > while the lighter fluid (density 
pi ^constant) is in the region y < 0; gravity is taken to point along negative y axis. 
As shown in Fig. 1, the magnetic field is taken along the z direction: 



II. BASIC EQUATIONS 



B = zBh{x,y,t); 



y>0 



(1) 



zBi{x,y,t); 



y<0 



so that V.B = 



(2) 
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automatically. The mushroom shaped perturbation interface which is called a bubble or a spike 
according as the lighter fluid pushes into the heavier fluid or the opposite is taken to have a parabolic 
form : 

y{x,t)=r]o{t) + r]2{t)x^ (3) 
Thus we have 

for a bubble: Vo > ^ ^72 < (4) 

for a spike: < and ri2 > (5) 

For uniform density fluid the equation of continuity V.v = is satisfled for irrotational fluid 
motion. Following Goncharov ^ the velocity potentials describing the irrotational motion for the 
heavier and lighter fluids are assumed to be given by 



y, t) = ai it) cos (A:x)e-'=(^-''» . ^ > q 



(6) 



U^.y.t) = 6o(t)2/ + 6i(t)cos(A;x)e^(^-''"W); y < 



(7) 



with Vh{i) = -V(f)h(i). 



(8) 



The fluid motion is governed by the ideal magneto hydrodynamic equations 



Ph{l) 



dv. 



dt 



+ (4{o-V)4{o 



-^Ph{i) - Ph{i)9 H (V X Bh{i^) X Bh(i) 

PHD 



(9) 



^ = V X [4(0 X 4(o] (10) 
For magnetic field of tlie form given by Eq. (1) 
^ (V X BhU)) X 4(0 = ^(4(0-V)4(0 - :T^^{B\ii)) (11) 



/^/l(0 /^/l(0 ^/^/l(0 

Substitution for ■u/i(/) from Eq. (8) in Eq. (9) followed by use of Eq. (11) leads to Bernoulli's 
equation for the MHD fluid 

^^Hi) , l^f7^ ^2 _ Phil) 1 p2 , A(0(^) ^9^ 

III. KINEMATIC AL AND DYNAMICAL BOUNDARY CON- 
DITIONS 

The kinematical boundary conditions satisfied by the interfacial surface y — rj{x,t) are 

| + W.| = W. (13) 

i'"h)xQ^ - = {Vh)y - {Vi)y (14) 

Prom BernouUi's Eq. (12) for the heavier and lighter fluids one obtains the following equation 

+i^-^)] + Mt)-Mt) (15) 



8 



— » 

Further with the help of Eqs. (1) and (2) and the incompressibihty condition V.fft^/) = 0, Eq. 
(11) simphfies to 

\- [Vh{i)-y)Bh{i) - U [lb) 

The interfacial kinematic boundary conditions (13) and (14) together with Bernoulh's Eq. (15) 
and magnetic induction Eq. (16) are employed in the next section to obtain the temporal evolution of 
the elevation of the tip of bubble (spike) like structures at the two fluid interface from its undisturbed 
level. 

IV. EQUATION FOR RAYLEIGH - TAYLOR INSTABIL- 
ITY INDUCED INTERFACIAL STRUCTURE PARAME- 
TERS 

Substituting ri{x,t) and (j)fi{i){x,y,t) from Eqs.(3),(6)-(8) in Eqs. (13) and (14) and expanding in 
powers of the transverse coordinate x up to i=2 and neglecting terms 0{x^){i > 3), we obtain the 
following equations 



= + 1)6 (18) 
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bo = - .J^' (19) 

- 2) 

bl = 7^7 TT«l (20) 

6 = ^?7o; 6 = ?72/A;; 6 = ^^oi (21) 

6 and 6 are respectively the nondimensionalized (with respect to the wave length) displacement 
and curvature of the tip of the bubble (spike) and ^3/k is tip velocity. 

At this stage it is in order to justify neglect of contribution from terms of order x\i > 3) as done 
here. This is provided on two counts: 

(i) The interface displacement y{x,t) is expanded in Eq. (3) keeping only terms of order a;^,-the 
customary practice in Layzer's approach. Since we are interested only in the motion close to the 
tip of the bubble or spike,i.e., for x it is sufficient to retain terms up to order x'^ and neglect 
0{x'){i > 3). 

(ii) Even if 77(2;, t) is expanded as 

r](x, t) = r]o{t) + r]2(t)x'^ + r]i{t)x^ + r]6{t)x^.... 

it can be shown that at the saturation level {drji/dt — 0) contributions from terms containing 
774, 776. ..are much smaller than that from r]2{t){see Appendix). Thus expansion of the kinematic con- 
dition and in its turn the expansion in Bernoulli's equation and Faraday's equation (which follows 
later) retaining higher order terms O(x^) can also be neglected. 
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Next let us turn to the magnetic field induction Eq. (16). To satisfy Eq. (16) with Vh given by 
Eq. (8) we set 

Bh{x, y, t) = Pho{t) + f3h{t) cos (A;a:)e-*^(^-''°W). y > q (22) 
in Eq. (16); this leads to 

Mt) + mt) + (3h{t)kf]o) cos (fca;)e-'=(^-^°(*)) - A;2ai/3^e-2'=(^-''°(*)) = (23) 

Corresponding to the parabolic interfacial structure represented by y{x,t) = ?7o(^) + V2{t)x'^ the 
foregoing equation yields on equating coefficients of {i — 0, 2) and neglecting terms 0(x')with 
i > 3 the following relation 

so that 

Pho{t) + f^h{t) = constant = Bho, say (24) 
i = 2: 

-^'~^,s^ SB.it) ^^-^ (25) 



5Bh{t) = 5Bh{t = 0) exp 



(6 + 2) . 



(26) 

so that 5Bh{t = 0) > (< 0); according as 5Bh(t = 0) > (< 0). 

In obtaining Eqs. (24) and (25) we have used the relation ^3 = ^1 = krjo = k'^ai (Eq. (17)). 

Similarly, to satisfy the magnetic field induction equation in the region y < 0, i.e., in the region 
occupied by the lighter fluid we set 

Bi{x, y, t) = Ao(t) + m) cos (A;a;)e'=(^-''°W); (27) 
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and proceeding as in case of the magnetic field induction Bh{x,y, t)m region y > 0, we obtain 
Ao(i) + Pi{t) = constant = Bio, say (28) 
and 



SBiit) Bio 

by using Eqs. (19) and (20) (^ bo + kbi + kai = 0). 

Again proceeding as in the deduction of Eq. (26) we obtain 



SBi{t) = SBi{t = 0)exp 



(30) 



/o''(6-|) (6-1) 

so that 6Bi{t = 0) > (< 0); according as 6Bi{t ^ 0) > {< 0). 

The magnetic field affected Rayleigh - Taylor instability induced growth of the mushroom shaped 
surface structure are determined by the parameters ii{t) , ^2{t) , ^3{t) as also the magnetic induction 
perturbation 5Bh{t) and 5Bi{t). To determine the time evolution of these five functions we need 
aside from the differential Eqs. (17), (18), (25) and (29) an extra one to complete the set. This is 
provided by Eq. (15). Now using Eqs. (22) and (24) one obtains 

^ Blix, y,t)^^- e^5Bnm2 + hx' (31) 



Similarly using Eqs. (27) and (28) one obtains 

:^Bf{x,y,t) = ^ + e^5Bim2 - hx' (32) 

Z/X; ZIJ,i III Z 

where 

\6Bh{t) \,\SBi{t) |<1 (33) 
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whenever the initial values | SBfi{0) \ and | SBi{0) |<S 1 as may be seen from Eqs. (26) and (30). 
This anticipation is substantiated later by numerical computation (Fig. 2 -Fig. 5). 

We next substitute for Bl{x,y,t)/2^h — Bi{x,y,t)/2jj,i from Eqs. (31) and (32) in Eq. (16) and 
use the dynamical boundary condition expressing balance of fluid and finite order magnetic pressure 
on two sides of the interface: 

Eq. (15) now reduces to 

P.h^ + ^(V0.)^] - Pi[-^ + ^(V0,)^] = -g{p, - p{)y + k^^5B,m, + i)a;^ 

+^dBim2 - bx' + hit) - flit) (35) 

which involves the influence only of the infinitesimal magnetic field fiuctuation on the interfacial 
structure. After some lengthy but straightforward algebraic manipulation we arrive at the required 
equation which is the last Eq. of the following set of Eqs. (36). The last Eq. of the set of Eqs. 
(36) is the required one as mentioned before and represents the dynamical boundary condition and 
obtained by setting y = rio + ri2X^ and equating coefficient of on both sides. All the equations are 
collected together below for the sake of convenience. 
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c I rr 



+ + . /7- (36) 



SB^it) (6-1) (6-1)^'/^ ^ 



dT D(6,r) (66-1) ' ' D{^2,r) 

(66 - 1) [,M55,(t)(2e, + 1) + ^55Kt)(26 - 1)] 



^(6,?^) ^ 9 



where,T = tJkg; r = ^; 1^(6, r) = 12(1 - r)e2' + 4(1 - r)6 + (r + 1); 
^ Pi 

iV(6, r) = 36(1 - r)e2 + 12(4 + r)6 + (7 - r) (37) 

V^(0 = \/Blo(^iQ)/ph{i)Ph{i) (38) 

is the Alfven velocity in the heavier (hghter) fluid. 

The above set of Eqs. describe the time evolution of a bubble. The time evolution of a spike 
is obtained from the same set by making the transformation 6 — ^ —6)6 — ^ —6 

r ^ lis. It 

is important to note that in the last Eq. of the set of Eqs. (36) the contribution to the bubble tip 
velocity ^ from the force of buoyancy g{ph — Pi)y is proportional to kg{ph — pi)^2 while that from the 
magnetic pressure fluctuation are proportional to k'^^Bf^/ iii)5Bi{t){^2 — \) and k'^^B'^Q/ ph)SBh{t){^2 + 
1) respectively as may be seen from Eqs. (31) and(32). Further both for bubbles and spikes 6 li^s 
in ( — 1, i); hence we always have (6 — ^) < and (6 + |) > 0. So by applying condition Eq. (30) 
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we find 

e^5Bim2-b<or>0; 

according as SBi{t = 0) > or <0. 

Similarly on applying Eq. (26) it follows that 

e^SBnm, + h>ar<0; 

according as 5Bh{t = 0) > or < 0. 



(39) 



(40) 



V. RESULTS AND DISCUSSIONS 

Analytical closed form solution of the set of Eqs. (36) not being feasible we take recourse to the 
method of numerical solution (5th order Runge-Kutta-Fehlberg method) and consider the following 
cases. 

Case A 

Assume Bho = 0, Biq ^ 0. Such a situation may occur when the lighter fluid (occupying the 
lower region y < 0) is ionized while the heavier fluid (region y > 0) is nonmagnetic. Prom Eqs. (32) 
and (36) and the concluding discussions of the foregoing section it is clearly seen that the instability 
driving pressure difference gi^ph — pi)^i is lowered or enhanced by | k -^SBi{t){$,2 — |)| according as 
5Bi{t = 0) is > or < 0. The concomitant growth rate modiflcations are shown in Fig. 2 (Fig. 3) 
which plots the bubble (spike) tip elevation | | and growth rate | I • Fig- 2 and Fig. 3 show that 
whether in case of suppression or enhancement the growth rate ^3(=| |) approaches an asymptotic 
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value as T — > oo both for bubble and for spike. This happens as SBi{t) exhibits similar asymptotic 
behavior as one may see in Fig. 2 and Fig. 3. The following analytic expressions for {^sjasymp as 
T — > oo are obtained by setting d^s/dr — together with Bho = in the last Eq. of the set of Eqs. 
(36): 



asympi bubble 



I 2Akg 



1 _ A kV"^ 



A 
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[(6) 



asympi spike 



I 2Akg 



1 - 2i 



l + A^kV,'^ 



[5Bi{oo)]spike 



(41) 



(42) 



A ' g 

Here 5i?i(oo) denotes the asymptotic value. The growth rate increases (destabihzation) if 5Bi{0) < 
(hence 5Bi{oo) < 0),i-e., the magnetic field perturbation diminishes the pressure below the interface 
relative to that above. On the other hand {^sjasym decrease and asymptote to (zero) as kVj^ /g 
increases if 5S;(0) > and therefore 5Bi{oo) > 0,i.e., the pressure below the interface increases and 
tends to restore stability. 

Case B 

Assume B^o ^ but Bi^ — 0. This situation is the reverse of that in case A and may arise when 
the heavier fluid is ionized while the lighter one is non magnetic. The dynamical boundary condition 
shows that following the same line of arguments as in case A but with but S/q = we find 

that the instability driving force g{^ph — pi)^2 is now enhanced or reduced by | k'^^SBh{t){^2 + |)| 
according as SBh{t = 0) is > or < 0. This conclusion is supported by the difference in the height 
of the bubble (or spike) tip shown in Fig. 4 and Fig. 5. However note that 5Bh{t) — > as i — > oo. 
This has the consequence that the asymptotic value of the velocity of the tip of the bubble (in spike) 
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height {C3)asymp IS the same as in the absence of magnetic field. But as | |= ^3 the height of the 
tip of the bubble (or spike) maintains a constant difference. 

Case C Assume both fluids are conducting and magnetic fleld is non zero on either side. We 
have considered two cases 

(i) If = ^ 5.0 

(ii) r^V;2 = 5.0, = 10.0 

(iii) rjlf = ^ 1.2 

(iv) r^n' = IVi' = 1-4 

with 6Bh{t = 0) = SBi{t = 0) > in each case. The bubble tip elevation ^1 as well as its velocity 
^1 = ^3 oscillates as the magnetic pressure acts on both sides of the interfaces but in opposite direction 
and with opposite phase. The results are shown in Fig. 6. In (i) and also in (iii) and (iv) the growth 
rate ^3 = ^1 oscillates approximately symmetrically about ,^3 = as r^V^ = ^V^ equation for in ^3 
in set of Eqs. (36); in (ii) the asymmetry results from difference in the driving pressure difference on 
two sides. Moreover it is to be noted from Fig. 7 and Fig. 8 that the oscillation frequency increases 
with 'k' and also with Alfven velocity. Occurrence of such an oscillation were also concluded for 
RMI [in]>[2l] increase in frequency similar to our case; however such oscillation are harmonic as 
against the nonlinear oscillations in our case. 

VI. SUMMARY 

Finally we summarize the results : 
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The change in interfacial pressure difference due to magnetic field fiuctuation leads to enhance- 
ment or suppression of instability as stated below. 

(i) If Bfio — 0,i.e, there exists no magnetic field above the two fiuid interface (y > 0) but 
Bio 7^ the lowering of the magnetic field below the interface y — due to an initial perturbation 
6Bi{0) < 0(^ 6Bi{t) < (by Eq. (30)) according to the fourth Eq. of the set of Eqs. (36)) leads 
to depression of pressure on the side of the lighter fiuid with the result that the instability growth is 
enhanced (Fig. 2). On the other hand if the initial perturbation SBi{0) > 0(=^ SBi{t) > (by Eq. 
(30)) the pressure on the side of the lighter fiuid increases with resulting suppression of growth rate 
which asymptote to (zero) as (-y-) increases (Fig. 3). 

SBi{t) tends to a constant value asymptotically as t ^ oo; this enables us to obtain an analytic 
expression for the asymptotic growth rate {^sjasym both for bubble and spike as given by Eqs. (41) 
and (42) respectively. 

(ii) If Bio = but Bho 7^ 0, it is 5Bh{t) — >■ (asymptotically whether initial perturbation 
5Bfi{0) > or < (Fig. 4 and Fig. 5). This has the consequence that the asymptotic growth rate 
becomes the same as in the nonmagnetic case. 

(in) If both Bfio 7^ Oand Bio with 5Bfi{t — 0) = SBi{t — 0) > the magnetic pressure 
perturbation acts on both sides of the interface but in opposite direction with opposite phase. This 
has the consequence that the growth rate ,^3 = .^1 oscillates symmetrically about ^3 = if = 
which increases in amplitude and frequency as the Alfven velocity increases. However, if ^ 
the oscillation is asymmetrical about ^3 = 0. 
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Appendix: 

Let the surface displacement rj{x,t) in Layzer's model expanded retaining higher powers of x:- 

r]{x, t) = r)o{t) + r]2{t)x^ + r)4{t)x^ + r}Q{t)x^ 

The time dependence of the coefficient functions rjiit) one obtained by equating coefficient of 
= 0, 2, 4, 6...) in the expansion of the kinematical condition in powers of x:- 

dri{x,t) _ 34, dr](x,t) _ _d± 
dt dx dx dy 

where the velocity potential (f){x, y, t) = ai{t) cos {kx)e~^^'^~'^°^*^'> 
This gives coefficient of x"^: 
^ = -ka, [3kri2 + ^' 

Coefficient of x^: ^ = -kai [dkrj^ - p'^rg - \k^r]2 - ^ 

Coefficient of x': ^ = -ka^ [ikrje - 7k%rj, + + + + f 

giving saturation values {drji/dt = 0): 
kr]2 = -k'^/6; 
krii = -fcVlSO; 
kr]6 = -A;V2835.... 
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Figure 1: Geometry of the model 
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Figure 2: Variation of ^i, ^2,bubble growth rate ^3(= ^i) and 6Bi with r for Vj^ = 
[Eq. |36]. Initial values ,^1 = 0.1,^2 = —0.1,^3 = 0.1 and r = 1.5 
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Figure 3: Variation of ,^2, spike growth rate ^3(= ^1) and 5Bi with r for = 
0[Eq.36] (with transformation ,^1 — )• —(,1,^2 —(,2,1" I/'" in Eq j36|) . Initial values 
= -0.1,6 = 0.1,6 = -0.1 and r = 1.5 
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Figure 4: Variation of .^i, ^2, bubble growth rate ^3{= C,i) and 6Bh with r for V^^ 



0[Eq.36 ]. Initial values ^1 = 0.1,^2 = -0.1,6 = 0.1, and r = 1.5 
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Figure 5: Variation of ^i, .^2, spike growth rate ^3{= C,i) and 6Bh with r for V^^ = 
OEq.36 (with transformation ^1 — )■ — .^,.^2 —^2,r — )■ 1/r in in Eq j36|) . Initial values 
^1 = -0.1,6 = 0.1,6 = -0.1 and r = 1.5 
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Figure 6: Growth rate oscillations for bubble. Initial values .^i = 0.1,^2 = —0.1, .^3 = 
0.1,6Bi{0) = 5BhiO) = 0.1 and r = 1.5 
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Figure 7: Oscillation of ^i, ^2, bubble growth rate ^3(= ^1), 6Bh and 6B1 with r as 
obtained by the solution of Eq. fl36|) with initial values ^1 = 0.1, ,^2 = —0.1, ,^3 = 
0.1,(55/(0) = SBhiO) = 0.1 and r = 1.5 



28 




10 20 30 40 50 60 



X 



Figure 8: Growth rate oscillation for bubble as obtained by the solution of Eq.f l36p 
with initial values = 0.1,^2 = -0.1,6 = 0.1,6Bii{0) = 55^(0) = 0.1 and r = 1.5 
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